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Certain approximate solutions of the continuum Schwinger-Dyson Equations (SDEs) predict chi- 
ral symmetry breaking in the 3d Thirring model when the number of fermion flavors Nf < 4.32 
^ whereas others predict symmetry breaking for all Nf [|[|. Our results from Monte Carlo sim- 
ulations with Nf = 6, predict a second order chiral phase transition. The critical coupling in this 
case corresponds to an ultra-violet fixed point of the renormalization group defining a non-trivial 
continuum limit. Further, our numerical simulations provide an estimate for the critical number 
of fermion flavors, Nf c ~ 6.5. 
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1. Introduction 

The study of quantum field theories in which the ground state shows a sensitivity to the number 
of fermion flavors Nf is intrinsically interesting. According to certain approximate solutions of 
continuum SDEs the 3d Thirring model as well as QED3 display this phenomenon. The Thirring 
model is a theory of fermions interacting via a current contact interaction: 

^ = +m ) Vi+ ^-(ifr iY ^ Vi )\ (1.1) 

where % are four-component spinors, m is a parity-conserving bare mass, and the index i runs 
over Nf distinct fermion flavors. In the chiral limit the Lagrangian of the model has a continu- 
ous U(l) chiral symmetry. Since the coupling g 2 has mass dimension —1, naive power-counting 
suggests that the model is non-renormalisable. However, as has been suspected for many years 
|], |5|], an expansion in powers of 1 /Nf, rather than g 2 , is exactly renormalisable and the model 
has a well-defined continuum limit corresponding to a UV-stable fixed point of the renormalization 
group. After the introduction of an auxiliary vector field, the above Lagrangian can be rewritten as: 



Sfan = Vi($+i£+m)y i + ^{A il ) 2 . (1.2) 

2g z ' 



The 1 /Nf expansion may not, however, describe the true behavior of the model, particularly for 
small Nf. Chiral symmetry breaking is forbidden at all orders in 1 /Nf, and yet may be predicted by 
a self-consistent approach such as SDEs ||]. SDEs solved exactly in the limit g 2 — ► 00 ^ show 
that chiral symmetry is broken for Nf < Nf c ~ 4.32. This is close to certain predictions of Nf c 
for the non-trivial IR behavior in QED3 Based on these results, at Nf c the model is expected 
to undergo an infinite order or conformal phase transition, originally discussed by Miranskii et al 
in the context of quenched QED4 [^]. Using a different sequence of truncations Hong and Park 
[§] found chiral symmetry breaking for all Nf. Previous Monte Carlo simulations [||, |9[ [l(| [jj]] 
provided evidence that for Nf < 5 the model is in the chirally broken phase at strong enough 
coupling. Here, we present preliminary results based on numerical simulations with Nf ■ = 2, 18 
in an effort to estimate Nf c . 

2. Lattice Model 

The lattice action we have used is based on [j|] and is as follows: 

S = ^L^W T 7M( x )( 1 + J ' A MW)^(^ + A)+h.c. 

xi o x\i 

ee Y,Xi(xMA,m}(x,y)Xi(y) + ^2ZM x ) 2 > (2-1) 

xyi ^8 x\l 

where %,% are staggered fermion fields and the flavor index i runs from 1 to N. We have intro- 
duced M[A,m] for the fermionic bilinear, which depends on both the auxiliary field defined on the 
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Figure 2: Chiral condensate vs. j3 for Nf = 6 near j3i; m . 



lattice links and the bare mass m. In three dimensions N staggered fermion describe Nf = 2N 
four-component continuum species. We employed the Hybrid Monte Carlo (HMC) algorithm to 
simulate even Nf and we used the Hybrid Molecular Dynamics (HMD) algorithm for odd or non- 
integer Nf. In the HMD simulations we used small enough fictitious time step At < 0.0025 and 
ensured that the 0(N 2 Az 2 ) systematic errors in the molecular dynamics steps were smaller than the 
statistical errors of the various observables. 



3. Critical Number of Fermion Flavors 

In this section, we summarize the discussion presented in regarding the non-conservation 



of the vector current at strong coupling in the lattice Thirring model of eq. ( 2. 1 ). At leading order in 
1 /N the vector propagator receives an extra non-transverse contribution from vacuum polarization, 
essentially due to the absence of the second diagram in Fig. [T[ It should be noted that in non- 
compact QED3 both diagrams are present and a cancellation of 0(a _1 ) contributions from the two 
diagrams in Fig. [j] occurs. In the 3d Thirring model however, the effect of the absence of the second 
diagram can be absorbed into a redefinition of the coupling: 

1 — g z J{m) 

where J(m) is the value of the integral contributed from the second diagram in the above figure. 
The physics described by continuum l/Nf perturbation theory occurs for the range of couplings 
g\ € [0,oo), ie. for g 2 <G [0,g 2 im ); to leading order in l/N l/g\ m = 2/3 for m = 0. In summary, 
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the strong coupling limit l/g 2 — > of the continuum model is recovered in the lattice model at 
\/g 2 — > g\ m > 0. As we already discussed in the Introduction, chiral symmetry breaking is absent 
in large-N calculations. Therefore, it may well be that the value of the second diagram in Fig. |T] 
is considerably altered in a chirally broken vacuum. Later in this section we will use this strong 
coupling lattice discretization effect as a criterion to decide whether chiral symmetry is broken or 
not at infinite coupling for different values of Nf. The data presented in Fig. || show clearly that 
the value of j8i; m (which corresponds to the infinite coupling limit in the continuum), where the 
condensate has its maximum value does not depend on the lattice size, because as mentioned in 
an earlier paragraph this is an 0{a~ l ) discretization effect. In addition, the same figure shows that 
within our statistical errors the position of the peak does not depend on the lattice size or the bare 
mass. Figure || shows that as expected the value of the chiral condensate (for m = 0.01) decreases 
as Nf is increased, because the fermion-antifermion screening has a stronger contribution for larger 
Nf. Furthermore, the value of j3i; m shifts to the left as the number of fermion flavors increases from 
Nf = 2 to Nf = 6. Figure |] shows that this trend is not monotonic; near Nf ~ 6.5, the curve of 
Ptim(Nf) reaches a minimum, implying a significant change in the strong coupling behavior of the 
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model. The minimum in py im (Nf) is a significant change in the phase diagram of the 3d Thirring 
model, implying that Nf c « 6.5. This conclusion is enhanced by the results shown in Fig. ^, where 
we plot the chiral condensate at pn m (Nf) versus Nf. The relatively rapid change in the value of the 
condensate around Nf = 6.5 may imply that a phase transition separates a chirally broken from a 
chirally symmetric phase with Nf c « 6.5. We are currently extending our simulations in order to 
understand better the properies of the chiral symmetry restoration transition. 



4. Analysis of Nf = 6 data 

During the last decade, Monte Carlo simulations provided evidence that for Nf < 5 and for 



sufficiently strong coupling the 3d Thirring model is in the chirally broken phase [|8j, g, |10|, 
In this section we present results from new simulations with Nf = 6 on lattice sizes varying from 
12 3 to 32 3 . As expected the phase transition for Nf = 6 is steeper than the transitions with smaller 
Nf. However, a detailed finite size scaling analysis of the Binder's Cumulant (not shown here) 
showed that the data at m = 0.01 are consistent with a crossover instead of a first order transition. 
In order to understand the properties of the model near the transition we fitted the data for the chiral 
condensate extracted from simulations on the largest 32 3 lattices to the Equation of State (EoS): 

m =A(r3-r3 c )(xx) p +B(xx) s ; (4.1) 

where p = 8 — 1/ /3 m . The assumption in eq. ( |4.1[ ) that finite size effects are negligible is justified 
by the observation that for the smallest mass used, m = 0.01, the difference between the values of 
(^W) extracted from simulations on 24 3 and 32 3 lattices is small. The fit was performed for /3 = 
0.34, ...,0.44 (see fig. ||), i.e. close enough to the transition but still away from j8/ im « 0.30 where 
strong coupling lattice discretization effects are severe. The results are: A = 1.38(13), B = 58(25), 
/3 C = 0.3082(36), 8 = 4.65(27), and j8 m = 0.28(2) with an acceptable fit quality # 2 / dof = 1.6. For 
Nf = 6 the value j8 f is close to j8/ !m , implying that Nf = 6 is close to Nf c , which is in accordance 
with the results presented in section 4. We also performed fits which included data from smaller 
volumes to a finite size scaling form of the EoS : 

mo =A((/J +CL-±)(XX) P + B{XX) 5 (4-2) 



5 



Three dimensional Thirring model 



Costas Strouthos 



0.2 r 



0.15 



Power Law fit to 32 

0.32<P<0.44 



0.1 - 



0.05 



o 


L=12m=0.01 





L=16m=0.01 


o 


L=24 m=0.01 


o 


L=32 m=0.01 


o 


L=32 m=0.02 




L=32 m=0.03 



n o 

" Si 
a S, 

o „ a 



°<>o -o 



0.1 



0.2 



0.3 
P 



0.4 



0.5 



Figure 6: Chiral condensate vs. j3 for Nf — 6. The dashed lines represent fits of the 32 3 , 0.32 < /3 < 0.44 
data to eq. 4. 1 . 



The results extracted from fits to eq. < ]4-2| ) (see fig. f7|) are of lower quality (% 2 /dof = 5.8) be- 
cause on smaller volumes the peaks near j8/, m are broader, implying that strong coupling lattice dis- 
cretization effects are more severe. In this case we got A = 1.91(4), B = 231(53), j8 = 0.3185(10), 
5 = 5.49(13), and j8 m = 0.19(2). 

5. Conclusions 

Our main results is that the Nf = 6 3d Thirring model like its smaller Nf counterparts appears 
to undergo a second order chiral symmetry breaking phase transition at strong coupling. Moreover, 
the Nf = 6 critical exponents are distinct from those of Nf = 2 and = 4 i.e. the different models 
have interacting continuum limits which are qualitatively, but not quantitatively similar. Further, 
a detailed analysis of the strong coupling behavior of the model for Nf = 2, 18 revealed that 
the critical number of fermion flavors where chiral symmetry is restored is Nf c sa 6.5 which is 
distinctly larger than the Nf c = 4.32 predicted by certain SDE approaches [[[]] and it contradicts the 
conclusion drawn from other SDE approaches which predict chiral symmetry breaking for all 
values of Nf. Our results may have broader interesting implications on the applicability of SDEs in 
non-perturbative field theories. There are predictions ^[|5]] that in the strong coupling limit the 
vector meson becomes massless suggesting that the IR limit of QED3 coincides with the UV limit 
of the 3d Thirring model. The problem of detecting dynamical symmetry breaking in lattice QED3 
is much harder than in 3d Thirring, because in QED3 the dynamical fermion mass is much smaller 



than the momentum cut-off [12]. 
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Figure 7: Chiral condensate vs. j3 for Nf = 6. The dashed lines represent fits of data from all lattice sizes 
for J3 > 0.34 to eq. 4.2. 
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